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OPTIMUM INTERCEPTION OF A BALLISTIC MISSILE 
AT INTiiiMii^IriTE RANGES 

by 

Frank D„ Faulkner^ 

Introduction o The problem studied in this paper is that of 
obtaining optimum trajectories, such as the one for effecting 
interception of a target following a known course above the sen- 
sible atmosphere, with minimum fuel consumption, or in minimum 
time« A procedure is given for determining the trajectory on a 
digital computer, and some conditions are given for the corres- 
ponding optimum trajectories. The methods apply generally to 
problems wherein the range and time are not too great. and the 
terminal velocity is irrelevant. 

An essential feature is in the use of the adjoint system of 
differential equations as defined by G, A, Bliss, A somewhat 
novel feature is in the use of direct methods and an optimizing 
principle in determining the trajectories, A simplification of 
the differential equations is effected by using the acceleration 
rather than the mass as a variable. 

The steering equation, which is well known, is a simple con 
sequence of the optimum principle in control theory which has 
received considerable attention in this country and in Russia re 
cently. The throttling relation, also obtained from an optimiam 
principle, is apparently not generally known: in problems where 

the energy available is limited, the energy input is to occur 
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when a vector defined by the adjoint system of differential 
equations is a maximm in magnitude. These principles reduce 
the problem to the solution of a three-point boundary problem. 

The following simplifying assumptions are made. The action 
takes place above the atmosphere so that aerodynamic forces are 
negligible. The trajectory of the missile is assumedto lie in 
a plane containing the center of the earth so that the trajectory 
is two dimensional, to reduce the number of variables. 



1. Basic equations. The equation of motion of a rocket in a 
gravitational field, subject to no outside forces, may be writ- 
ten as 

[1] ? = g + a, 

where r is the position vector, g is the acceleration due to 

gravity, a is the acceleration due to thrust, and a dot (°) ova? 

a variable indicates its time derivative. We may write 

[2] a = c^me/Cl-m), 

where m is the ratio of the mass of fuel which has been consumed 
to the initial gross mass of the rocket e is the unit vector 

in the direction of the thrust^ and c ' is a constant, 

A useful kinematic relation connecting the acceleration and 
the fuel consumption is the following 



[3] 




-c'ln(l-m) , 



Since the fuel consumed is prop^tlonal to m and In(l-m) is a 
mono tonic function of m, conditions involving the final value of 
the mass may be rephrased in terms of the integral of the accelera- 
tion, subject to the constraints on the size of m. For practical 
purposes we may take m as bounded above by some constant m 
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and bo\onded below by zero, and it may be chosen anywhere on or 
between these limits, so long as any fuel remains. 

The equations of motion then have the form 
,x = gi + a cos p 
y - g2 + a sin p , 

where we may think of gj. as -c x/r® , with c a constant and 

o o 

p 2 = 35 - 2 + y2 j and gg as -c y/r®; these equations may be refined 

o 

or simplified. 

To get the needed formulas, let us multiply equations [4] 
through by two new variables u,v, which are unspecified so 
far, and integrate formally to get 
„T 

[5] [u(x -gi -a cos p) + v(y -gg »a sin p) ]dt = 0 , 

'V 

where x,y are any solutions to [4], 

Let us consider also a neighboring path whereon the original 
equations are satisfied. For the first variations, equation [5] 
becomes 



T 

[6] I [u(6x - gi^fix »gi„6y » 6a cos p + a 6p sin p) 

Jq X y 

+ v(6y - g 2 x^x »g2y5y - 6a sin p - a 6p cos p) ]dt = 0; 

subscripts denote partial derivatives g^^ = 0gi/Sx, etc. If 
we integrate by parts, this may be rewritten 



[7] [u6x - u6x + v6y - v6y!^+ [a(u cos p + v sin p) ] 
T 

+ J^C^X(U -gi^U “g2x^) + " glyU - ggyV) 




1 



- 6a(u cos p + V sin p) 

- a 6p(-u sin p + v cos p) ]dt = 0; 



here t^ is a symbol for any time or times when a is 
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discontinuous. 

Now, to simplify equation [7], let us choose UjV as solu- 
tions to the system of differential equations 



this choice of u,v eliminates the variations of the dependent 
variables x,y from the integral in [7], Equations [8] are 
called adjoint to the variational equations corresponding to 
[4]; in general, it is an integration by parts and setting to 
zero the coefficients of terms involving the dependent varia- 
bles which determines the adjoint system. Let us further choose 
the solutions so that u(T) = v(T) = 0, Equation [7] is then 
a fairly general relation connecting the variations of the end 
values of x,y with the variations of the control variations 
a,p for any particular fixed value T and any particular trajec- 
tory. 

2, Variations of dependent variables. If we further choose the 
solution Ui,Vi to [8] such that 



and if all of the initial values are given, equation [7] becomes 




oo _ 

-u “ - g2 V = 0 



o 



[9] 




0, 




+ [6a(uj,cos p + Vj^sin p) + a 6p(-UiSin p + v.cos p) ]dt 
*^0 



similarly, if U 2 (T) = 0, V 2 (T) = -1, 

[11] 6y(T) = -[a(u 2 Cos p + V 2 sin + 
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+ [6a(u2Cos p + VgSin p) + a 6p(-U2sin p + V2C0S p)]dt, 

t 

These are the essential equations for variation and control, 
since they express the effect on the terminal values of x,y 
of small changes in the control variables a,p. 



It will be necessary to have a fundamental set of solutions 

uationj 
3 .3 4 4 



112 2 

to equations [8]; we may choose, to be definite, u ,v ,u ,v , 



u ,v ,u ,v , such that 



[ 12 ] 



.u^(0) = 1, v^(0) = 0, u^(0) = 0, v^(0) = 0 
'u^(O) = 0, v^(0) = 1, u^(0) = 0, v^(0) = 0 



u (0) = 0 , ® " " 

^u"^(0) = 0 ,»->» v^(0) = 1 . 

Then, for any particular trajectory and any time T, every 
solution is a linear combination of these, and we may write 

/“i = i =13“^ 

[13] J I 



i = 1,2 



for the proper choice of the c’s< 



3o Maximizing principle. It is a property of extremals in nor- 
mal problems that they furnish a maximum to an integral and that 
this is done by maximizing the integrand; this may be taken as 
the characteristic of an extremal as follows. 

Let us write [5] in the alternate form 

[14] [ux - ux + vy - vy]J? + [ (3^ - ug^ + yv - vg2)dt 

u .1^ 




f 
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Now suppose we have selected w = ui + vj in some way, where 
u,v satisfy [8] and w = |w(t)| is a diminishing function of 
t with w(T) = 0« Without regard for a reason, let us consider 
the problem of maximizing the integral on the right in [14], by 
choosing a properlyo It most satisfy the constraints that (i) 
0 ^ a ®taax’ ^max ^ known function of t, and (ii) 

pT 

[15] a dt s C, a given constant; 

^0 



C is by [3] a measure of the fuel consumption. For a^^^ we may 
take the function 

[16] a ov = l6sser[c'm__ /(l“in„„^t) , cm /(1-m )], 

'• max ’• max: max ’ max^ max ’ 

where ^max’^max maxinumi values of m,m, determined 

by the rocket mechanism. 

We see that to maximize the integral in [14] it is neces- 
sary and sufficient that (i) a is parallel to w, so that 



[17] 



tan p = v/u 



(the proper branch of p chosen) and (ii) 



[18] 



a = 



max 



0 ^ t ti 



0 , ti ^ t s T, 



ti being chosen so that equation [15] is satisfied. 

An outline of the proof is given here; it is given in 
detail in (1, section 5), Consider the path defined above and 
any other path; let the functions for a,p be denoted by A,P 
on the second path, A must satisfy the constraints (i),(ii) 
specified on a above, so that 



[19] 



f A dt ^ C< 
‘’O 



It is seen that 




- J A w cos(p-P) dt 0 
■^1 

The s\am of the first and third terms is positive or zero and 
the second term is positive or zeroo This establishes that the 
path just described maximizes the integral, 

4. Attaining a fixed point in specified time with minimum fuel. 
Let us consider, as an example, the problem of determining a 
trajectory to attain a specified point X^,Y^ in a specified 
time T with a minimum of fuel, assuming that all of the initial 
conditions are given. We will first take up a routine for 
finding the trajectory and then show how it satisfies certain 
conditions which can be checked. 



where 0 is a number to be fo\md. Each of the quantities should 
have an iteration mmiber, since they will be in an iterative com~ 
putational routine. Now let 



properly chosen. Let us guess also tjj maximum thrust is to be 
applied until time t^, and thrust is zero thereafter. 

Now compute the corresponding trajectory. It will lead to 
terminal values x(T),y(T), which are in error by amounts 
X^-x(T) ,Y^-y(T) , Compute simultaneously a fundamental set of 
solutions u^,v^ to the adjoint system of differential equations 
[8], and two further integrals which occur later in equation [22], 



Let us guess initial values for the solutions ui,Vi and 
solutions to [8], and a linear combination of these 




[17] tan p = u/v 
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From [17],[20], for fixed or determined Uj. jVj jUg ,Vg , we get 



[21] 6p = 6 ©(UiV2 - U2Vi)/([u]2+[v]2) „ 

The relations for the variations of the end values then become 



[ 22 ] 



6x(T) 



V6y(T) 



= [aCuiCos p + v^sin p) 

p^i 

+ a(-uisln p +^V|_Q_QS^)(ui^^^ 

= [a(u2Cos p + V2Sin + J"”” 



Now set 

.6x(T) = X„ - x(T) 

[23] { ^ 

^6y(I) = Yf - y(T) , 

substitute into [22], and solve for 6ti,6©c. This gives corrected 
values of t^,© for the next round, 

1 4 

From the fundamental set of solutions u , ® » »v , corrected 
estimates are made for the initial values of Ui,Vi,U 2 ,V 2 and 
their derivatives. All of the initial values are then available 
for starting the next iteration. 

This routine is continued until some convergence criterion 
is satisfied; for example one might require 

[24] E = (X^ - x[T])2 +(Y^ c y[T])2 
to be less than some preassigned number. 

If the routine given above converges, we have found a 
trajectory which effects at least a stationary value for the 
fuel consumption compared with paths in some neighborhood. For 
if we consider neighboring paths in equation [14], we see that 
the first variation of the left side vanishes, by virtue of 
equations [8] and the fact that 6x(0) ,6y(0) ,6x(0) ,6y(0) ,6x(T) , 
6y(T) ,u(T) ,v(T) all vanish. But the right side of the equation 
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Is positive for every admissible variation of the acceleration, 
as was shown in section 3o Hence there is no neighboring path 
whose first variations do not vanish which requires the same 
or less fuel. If gravity can be approximated as a linear func- 
tion of position, this also establishes the sufficient conditions 
for an absolute minimum [see (1), section 6], 

5, Variations of terminal time; minimum fuel consumption. Now 
suppose that the terminal time T may vary. The change in the 
end values of x,y are 



,6X = x(T)6T + 6x(T) 

[25] { 

''6Y = y(T)6T + 6y(T) , 

There are corresponding equations for the terminal values of 
Ui, etc,, 

Ui(T) = -Ui(T)6T, Vi(T) = -Vi(T)6T, Ui(T) = -1 -u\(T)6T, 
etc,; in these equations T should have an iteration subscript, 
being the approximation associated with the last trajectory run. 
Since there are three variables ti,T,0, and only two equa- 
tions for X^,Y^, a further relation is required. If we are in- 
terested in minimum fuel consumption, then t^ must be a minimum. 
Equation [7] may be rewritten 

[26] (-u6x - v6y)j = [a(u cos p + v sin p)]|. _6ti , 
if terms that must vanish are omitted. 

If we consider two neighboring trajectories, each of which 
effects interception of a target whose coordinates are X^(T), 
Y^(T), then 6x(T) = (X^-x)6T, 6y(T) = (Y^-y)6T, and [26] 
becomes 

[27] (“W) j»[(jt^-x)I + C^|.-y) j]j6T = (aw)^^_6ti. 








§«P 
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Hence, for minimum fuel consumption, 

[28] w(T) °[(X^-x)I + = 0: 

these two vectors must be perpendicular; apparently this is 
a fairly general condition associated with minimum fuel con- 
sumptiono This condition allows a direct determination of the 
angle Oo This in turn defines 6© as the difference of the 
two successive values of ©o 

The value for 6x(T) which must be used in the iterative 
routine is 6x(T) = X^(T) - x(T) + - x)^6T, and the equa- 

tions for corrections to T,tj are 

,X.(T) - x(T) = - (X. - x)6T + 6x(T) 

[29] { ^ , 

^Y^(T) - y(T) = - (Y^ - y)6T + 6y(T) , 

where 6x(T),6y(T), linear in 6©,6ti, are given by equations 
[22], Since 6© is already determined, these equations [29] 
then yield 6T,6tjo 

The iteration is continued until some convergence criterion 
is satisfied. It must also be checked that w is a decreasing 
function of time; actually all that is required is that 
w(t) > w(ti) when t < tj and w(t) < wCt^) when t > t^. It 
should be pointed out that the relation involving the end con- 
dition, equation [28], is only a necessary condition. There do 
not seem to be any simple ways to ensure that a minimum, rather 
than a maximxam is obtained; in the case of constant gravity and 
a ballistic target, the trajectory always furnishes a minimum, 

6, Comments, The procedure seems to furnish a rather simple ex- 
tension of the procedures which Bliss (3) introduced in Ballistics, 
for calculating differentials. There seems to be only one 
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peculiarity of the trajectories, namely that the vector w 
is decreasing in magnitude o 

If a minimum time of interception is desired, then set 
tj = lesser( o Equations [29] then furnish the changes 
6T,6©o It is necessary to make the first estimate of © in a 
certain range, else the routine may not converge, or may 

converge to the trajectory which takes maximum time (l,p 17) » 

The use of the acceleration rather than the mass as a varia- 
ble simplifies the differential equations, since the mass enters 
non linearly in the differential equations o On the other hand, 
the bound on a is a function of time, not a constant, in the 
above example, and in more general cases it does not seem possible 
to express the bounds explicitly as functions of time. 

There seem to be two general principles associated with 
optimimi trajectories where there is a limited amount of energy 
available, so that throttling must occur in some form or another „ 
The first of these is that the orientation of the energy-input 
vector, the acceleration in this case, is such as to maximize 
its projection onto the adjoint vector. The second is that the 
energy is to be expended at a time when the adjoint vector is 
large. The first condition above leads to the well-known steering 
criterion, apparently first published by Lawden [see (2) for refer 
ences]. The second is apparently equivalent to the Weierstrass 
condition in some problems, though as expressed above, it has no 
analog in classic calculus of variations. 

The method of solution, using procedures due to Bliss (3), (4) 
(6), has the nice feature that, should the actual trajectory 
deviate in some way from the planned trajectory, the same equa- 
tions for the variations can be used to determine a corrective 
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thrust schedule,, The vectors Wi,W 2 may be considered the 
impulse response vectors for x(T),y(T) respectively, since 
the dot product e<>w^ gives the change in the terminal value 
of the corresponding variable due to a unit impulse in the 
direction of e o 

In the computing which has been done, a modified Runge- 
Kutta routine was used which has the desirable feature that it 
is easy to change the time intervals, which must be done in the 
neighborhood of ti,T„ No convergence problems have been en- 
countered in the ballistic missile interception problemo On 
the other hand it is apparently necessary to make a close initial 
estimate in a corresponding orbital transfer problemo The values 
obtained from the case where g varies at most linearly should 
furnish reasonble starting values, should convergence be a problemo 

Considerable attention has been attracted recently by the 
maximum principle, due largely to papers by Pontriagin (6)o To 
the author this seems to be the fundamental way to approach prob- 
lems in the calculus of variations; the Euler equations 

are derived by effecting a stationary value for an integral 

through the e6 processes of calculus. In most problems there 
seems to be no practical new information in the maximum prin- 
ciple. Superficially the system of differential equations is 
of lower order when the maximum principle is taken as the basic 
concept. In the nmerical solution, the last Euler equation [17] 
may be solved by Newton's method, which depends on derivatives 
so that it is equivalent to an additional differential equation. 

Several authors have dealt- with the mass and handled the 
bounds on m by introducing a new variable, such as (j) in the 
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equation 



m 



m. 



max 




In simple problems, both the mathematical theory and the pro- 
gramming of the numerical solution are simpler in terms of the 
original variables; the logical decisions are of a type built 
into all general-purpose digital computers „ 

It seems to the author that one of the most important 
features in recent control theory rests on the use and inter- 
pretation of the adjoint system of differential equations o 
As remarked earlier, the fomdations were laid by Bliss, with 
more recent applications and contributions by Drenick (7) , 

Tsien (8), Breakwell (9), and Tyndall (10), among otherSo The 
mathematical contribution to calculus of variations is in the 
n\americal procedures for solving two- (and more) point boundary- 
value problems on high-speed computerso Another method of solu- 
tion, the gradient method of Kelley (11) also utilizes the 
adjoint systemo It is sometimes stated that the adjoint equa- 
tions [8] are necessary conditions for an extremalo This seems 
a logical error; the variables u,v are introduced by the 
mathematician, not implicit in the problem, and in turn are 
chosen to satisfy [8], because that choice simplifies the inte- 
gral in equation [7]o Dro So Ross usually brings out this fact* 
In problems such as the above, the construction of the functional 
of Pontriagin [(6), p 16] is equivalent to solving the three- 
point problem involving the original system of equations and 
the adjoint systemo 
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